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Abstract 

"-^ LuCuGa04 has magnetic Cu 2+ and diamagnetic Ga 3+ ions distributed on a triangular bilayer. 

Susceptibility, /iSR and neutron scattering measurements show that at low temperature the spins 
form a short range correlated state with a gapless continuum of excitations. The development of 

T-H 

this spin liquid state is strongly reminiscent of that in other two dimensional frustrated magnets 
£ — . and we suggest that it is a generic property of such systems, even surviving in the presence of 

strong quenched disorder. 
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I. INTRODUCTION 



LuCuGa0 4 encompasses several topics of interest in condensed matter physics. It has 
the same crystal structure as the multiferroic LuE^Oj^ both contain triangular bilayers 
(Fig. [TJ, a two-dimensional geometrically frustrated lattice which has so far not been widely 
studied. Theories of resonating valence bond (RVB) states and other spin liquid ground 
states 2 typically address low-dimensional frustrated magnets with s — 1/2 and experimental 
realizations are sought-after for comparison. In LuCuGa04 the magnetic Cu 2+ ions have 
s — 1/2, suggesting the possibility of quantum magnetic phenomena, but are strongly diluted 
by diamagnetic Ga 3+ giving a material with quenched disorder and charge frustration^. 

At a classical level, spins with antiferromagnetic couplings residing on lattices with trian- 
gular geometry are frustrated: all pairwise interactions cannot be minimized simultaneously. 
Many other systems map onto such models, for example the configurations of two cations 
on such a lattice are identical to those of antiferromagnetically coupled Ising spins on the 
same lattice^, a situation which is known as charge frustration. In LuCuGaO^ the Cu 2+ and 
Ga 3+ cations are distributed over a single site forming well separated triangular bilayers. A 
single bilayer contains two opposed triangular lattices (see Fig. [T| and the shortest distance 
between cations is between opposite faces of the bilayer, forming a puckered honeycomb lat- 
tice. The combination of honeycomb and triangular geometries may be conveniently viewed 
as a Ji — J2 model on the honeycomb lattice^, which is a frustrating geometry for magnetic 
order (see Fig. [TJ. LuCuGa0 4 is therefore expected to exhibit both charge and magnetic 
frustration. 

The combination of a frustrating lattice geometry and s = 1/2 is expected to lead to 
unconventional groundstates such as a spin liquid, as the energy can be further reduced 
below the classical groundstate by quantum fluctuations (many types and definitions of spin 
liquids are reviewed by Normand in Ref. [2]). Whilst geometrically frustrated magnets have 
typically been viewed as clean systems in which the effect of frustration can be isolated (in 
contrast to spin glasses, where structural disorder is also inevitable), it has become apparent 
that weak structural disorder or non-magnetic impurities can determine the low temperature 
behaviour of frustrated magnets, where any perturbation of a fragile groundstate is magnified 
in importance. Examples include the pinning of dimers around non-magnetic impurities in 
the s = 1/2 kagome Heisenberg antiferromagnetP^or the strain- induced spin glass transition 
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FIG. 1: The topology of the triangular bilayers formed by Cu + and Ga + cations in LuCuGa04. 
The two faces of the bilayer are shown in black and white. The shortest bonds (called J\ in the 
text) are those that link the two faces of the bilayer forming a puckered honeycomb lattice (grey). 
Although the honeycomb lattice is bipartite and cation order is therefore possible, the side views 
(middle and bottom) show that the puckering of this lattice means that such order would cause 
charge separation across the layer. The bond distance within the faces of the bilayer are slightly 
longer and form triangular lattices (black and dashes) on which cation (and magnetic) order will 
be highly frustrated. 

in the pyrochlore Heisenberg antiferromagnet with weak bond disorder^. 

LuCuGa04 is a member of the series RMM'O^, where R is a small rare earth such as Yb 3+ 
or Lu 3+ , M is a late 3d transition metal cation with +2 charge and M' is a +3 cation from the 
late transition metals or main group. Much attention has been focussed on LuFe204, where 
charge ordering of Fe 2+ and Fe 3+ is possible and is thought to underlie the multiferroicity^. 
Numerous examples with two different cations have also been characterized. As well as 
LuCuGa0 4 , YMnFe0 4 , YbMFe0 4 (M = Mg, Fe, Co, Cu), YbCuGa0 4 , LuMFe0 4 (M = Zn, 
Fe, Co, Cu) and LuCoGa0 4 have all been studiecP™ Structurally, all are characterized 
as described above, with randomly distributed M and M' ions on the bilayer (the crystal 
structure is illustrated in Fig. [2] and described in greater detail below). With the exception 



of i?CuGa04 (R = Yb, Lu)^ all these materials exhibit a splitting of field cooled and zero 
field cooled magnetic susceptibilities at temperatures of order 20 K (and in some cases 
considerably higher), indicative of freezing or spin glass transitions. Previously the most 
extensively investigated stoichiometry was LuMgFeO j 11 * 18 * ^^. Studies using magnetization 
measurements, neutron scattering and Mossbauer spectroscopy established that there are 
short range magnetic correlations associated with frustrated two-dimensional ordering and 
suggested that random dilution leads to the formation of clusters and two types of site: 
those lying in the body of clusters and having a full complement of magnetic neighbors, and 
those lying on extended branches of clusters and having just one or two neighbors. 

Previous studies of LuCuGa04 also suggest that the copper and gallium ions are randomly 
distributed on the sites of the bilayer and consequently that LuCuGa04 is both spin and 
charge frustrated^. Indeed, despite a Curie- Weiss temperature of 6cw = — 69 K, no magnetic 
ordering or freezing transition is observed down to 0.4 K where a broad peak in the ac- 
susceptibility is currently attributed to a spin glass transition^. The primary difference 
between those RMM'O^ with high temperature freezing transitions and LuCuGaC^ is the 
presence of (only) s = 1/2 moments. In this paper we use polarized neutron scattering, /iSR 
and neutron spectroscopy to produce a detailed microscopic characterization of LuCuGa04 
which suggests that such a state exists, even in the presence of strong structural disorder. 

II. EXPERIMENTAL METHODS 

A 20 gram powder sample of LuCuGa04 was prepared by reaction of the appropriate 
stoichiometric quantities of Lu 2 3 (Alfa Aesar 99.99%), CuO (Sigma Aldrich, 99.995%) and 
Ga 2 03 (Alfa Aesar 99.999%). CuO and LU2O3 were pre-treated by heating in air for 24 hours 
at 900 and 1200 °C respectively. The starting materials were ground together, pelletized 
and heated in air at 1050 °C for a total of 120 hours with 1 intermediate grinding. The field 
cooled and zero field cooled (FC/ZFC) dc-susceptibility was measured between 1.8 and 300 
K in a field of 1000 G using a Quantum Design MPMS-7 SQUID magnetometer. 

Neutron powder diffraction patterns were obtained at 70 K and 1.5 K using the D1A 
diffractometer at the ILL with an incident wavelength of 1.91 A and standard Orange cryo- 
stat. Polarized neutron scattering using the xyz technique was performed on DT^at the ILL 
using a wavelength of 4.8 A at temperatures of 0.08, 0.5, 5 and 50 K. Low temperatures were 
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FIG. 2: Crystal structure of LuCuGaC-4 (the c-axis is vertical on the page, the projection axis 
is the a — b diagonal) as obtained from the Rietveld refinement shown in Fig. [3| The bilayers 
(black bonds) are formed of edge-sharing M^Os trigonal bipyramids are shown in green, 

2 ~ in red). The bilayers are separated by layers of edge-sharing LuOg octahedra. The Lu 3+ 
cations randomly occupy one of two sites displaced from the centre of each octahedron (blue). We 
suggest (see discussion) that the minimization of the Coulomb interaction amongst the cations is 
frustrated, and that short range cation correlations are formed amongst the Cu 2+ and Ga 3+ on the 
bilayers. The siting of a particular Lu 3+ ion would also be controlled by the cation distributions 
in a triangle on the face of the bilayer immediately above and below the LuOg octahedra. 

achieved by using a dilution refrigerator insert in a standard Orange cryostat. The sample 
was contained in a copper can with 3 He exchange gas to provide thermalization below 1 K. 
An inelastic neutron scattering experiment was carried out on the direct geometry chopper 
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spectrometer IN4 at the ILL, at temperatures of 1.6, 5 and 50 K using an incident energy 
of 17 meV. 

A zero field and longitudinal field muon spin rotation study was conducted between 0.06 
and 50 K with fields up to 2500 G using the MuSR instrument at the ISIS pulsed muon 
source. The sample was pressed tightly in a silver disk-shaped holder approximately 2 mm 
deep and cooled with a dilution refrigerator. 

III. RESULTS 

A. Characterization 

We have used powder neutron diffraction and magnetic susceptibility measurements to 
establish the basic quality of the sample and compare to previous studies. The powder 
neutron diffraction data show that the sample is phase pure and highly crystalline. The 
Rietveld method was used to refine three structural models as described in Ref. IU (see also 
Ref. HI]). The crystal structure has the space group R3m and consists of double layers of 
edge-sharing M^'^0 5 triangular bipyramids connected by triangular layers of ROg octahedra. 
All three models have Cu 2+ and Ga 3+ randomly distributed on the 6c site, and differ with 
respect to the position of Lu 3+ . The R 3+ ions lie at the centre of distorted octahedra of 2 ~ 
ions, directly above and below the centre of triangles of transition metal ions. The simplest 
model has Lu 3+ positioned at (0, 0, 0), midway between adjacent bilayers. The more complex 
models allow anisotropic motion of the Lu 3+ ion, or a shift to a lower symmetry site with 
random half occupation at (0, 0, ±z) (z = 0.009 ±0.0002, corresponding to a displacement of 
0.22 A). In this case, the position of each Lu 3+ is thought to be controlled by the distribution 
of cations on the coordinating triangles of adjacent bilayers^. As in Ref. H]the third model 
provides the best combination of fit and realistic parameters. There is no indication of any 
structural change between 70 K and 1.5 K. The crystal structure is illustrated in Fig. |2j 

We have fitted the inverse magnetic susceptibility to the Curie- Weiss law between 50 K 
and 300 K, obtaining 8qw — —62.1 ± 1.8 K . This is somewhat smaller than the previously 
reported value (-69 K), but in Ref. Ulthe fitting range appears to be restricted to 50 K to 
150 K. The magnetic moment obtained from the Curie constant is 2.2 ± 0.03 fiB per Cu 2+ 
(the expected value for S = 1/2 is 1.7 /ib atom -1 ), and the effective moment at 295 K is 
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FIG. 3: Neutron powder diffraction at 1.5 K. The experimental data are the red crosses. The 
fitted line (green) was obtained by using the Rietveld method to refine the model M3 of Ref. HI 
The corresponding structure is shown in Fig. [2j The calculated peak positions are indicated by 
the black ticks and the difference between experimental and fitted profiles by the magenta line. 

fi eS = v^S^T = 1.73 ± 0.06 /xb per Cu 2+ . The effective moment falls at low temperature 
where the inverse susceptibility becomes non-linear, indicating increasing antiferromagnetic 
spin correlations. At these temperatures there is no discernible splitting between field cooled 
and zero field cooled measurements. 

We have also fitted the entire temperature dependence of the inverse susceptibility with 
a model of "orphan" spins^. This is essentially a Curie- Weiss law with an additional Curie 
contribution due to isolated spins, which is expected to become dominant at low temperature. 
Since these spins are, in principle, isolated, they should have 6 C w — 0. In this case Qcw\ — 
-67±2.9 K, m = 2.1±0.008 /i B , 9 CW 2 = 1.5±0.8 K and /i 2 = 0.33±0.06 /i B - The downward 
curving form of x -1 is a rather general feature of many frustrated antiferromagnets (see 
Ref. [28] and references therein) and the reality of the orphan spins remains debatabl e! 27 * 28 *. 
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FIG. 4: Susceptibility and (inset) inverse susceptibility for FC dc SQUID measurements at 1000 
G. The solid line is a fit to the Curie- Wiess law between 200 and 60 K. ZFC results are identical 
to FC within experimental error. Also shown is the fit to a model with "orphan" spins^. 

B. Polarized Neutron Scattering 

The XT/z-polarization analysis technique was used to separate nuclear, spin incoherent and 
magnetic contributions to the total scattering^. The temperatures were selected to compare 
the paramagnetic regime, where the susceptibility is still described by the Curie- Weiss law 
(50 K), with the regime in which the inverse susceptibility becomes non-linear (5 K), and 
the low temperature regime around the proposed spin glass transition. Furthermore, the 
separation of the contributions allows us to distinguish the origin of any short range order 
(i.e nuclear or magnetic). The nuclear channel (not shown) contains only Bragg scattering. 
There is a clear change in magnetic scattering between 50 K and all lower temperatures 
studied. However, there is no significant change in the low temperature scattering from 
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5 K down to 0.08 K, in particular between 0.5 and 0.08 K across the proposed spin glass 
transition, so we have combined the data sets from 0.08, 0.5 and 5 K. As shown in Fig. |5j 
at 50 K the magnetic scattering can be approximately fitted by the form factor of Cu 2+ , 
as expected for a paramagnet. At low temperature the form factor response at small \Q\ is 
suppressed and a broad feature appears, centered around 1.25 A -1 . There is some additional 
scattering at small \Q\, above the fitted form factor, which we attribute to air scattering. 
The magnetic scattering is very weak and has therefore not been observed using conventional 
neutron powder diffractometers (there is no sign of this feature in Fig. [3] or the earlier study 
of LuCuGaOj 1 ). 

Using an expression for the structure factor of disordered near neighbor antiferromagnetic 
dimers, 

a fit was obtained, as shown in Fig. [HJ with F 2 (Q) being the form factor for Cu 2+ and d the 
distance between the AFM correlated Cu 2+ . No appreciable difference in the fit is observed 
if the value of the Cu-Cu distance is 3.44 A (the distance between Cu-Cu (or Cu-Ga etc) in 
the triangular layers), 3.04 A (the Cu-Cu distance across the bilayer), or a weighted average 
of the two (i.e. due to neighbors on both faces of the bilayer). 

This form suggests that there are magnetic correlations between nearest neighbor Cu 2+ 
ions, which develop in the non-linear region of x~ X - The position and form agrees with other 
members of the series where the larger moments of Fe 3+ or Co 2+ render the diffuse magnetic 
scattering more readily observable^- 1 . In LuMgFe04, the dominant feature is also a peak 
at Q ~ 1.25 A -1 , but the diffuse scattering profile is much more detailed and correlations 
could be detected up to the twelfth neighbor shell. This fitting strategy involving further 
neighbor correlations^ produced unphysical values for further neighbor correlation functions 
in LuCuGaOj so we conclude that the correlations are very short ranged. 

We gain no particular understanding of the anomaly previously observed in the suscepti- 
bility and specific heat at 0.4 rd^as there is no difference between the data measured at 0.08 
K and 5 K. This is also a feature of LuMgFe04 where no change in the neutron scattering 
intensity is seen across the susceptibility cus p 118 1 21 1 . 
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FIG. 5: Separated magnetic component from xyz polarization analysis experiment. The combined 
low temperature results (0.08 to 5 K) fit well to a model of nearest neighbor antiferromagnetic 
correlations. At 50 K LuCuGaO"4 is paramagnetic and the scattering approximately follows the 
magnetic form factor for Cu 2+ ^. 

C. /xSR 



The development of the spin correlations at low temperature was further investigated 
using /iSR. The decay of the implanted muon polarization was measured in zero field as a 
function of temperature and as a function of magnetic field at selected temperatures. /iSR 
is sensitive to magnetic fluctuations with nuclear and electronic origins on a timescale of 
10 -5 - 10~ 7 seconds. In a paramagnet the electronic fluctuations are mostly very rapid and 
lie outside the window of the /zSR experiment. Normally in such a situation, any observed 
relaxation is a consequence of the nuclear component. The nuclear component of the muon 
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relaxation can usually be decoupled with a field of the order of 10 — 100 G. We have applied 
longitudinal fields at several temperatures. At 50 K, where the polarized neutron scattering 
also shows the sample to be paramagnetic, the /iSR signal is decoupled in fields of less than 
100 G, suggesting the signal is entirely due to relaxation of the nuclear spins. However, at 
low temperatures, where the polarized neutron scattering shows that spin correlations are 
developing, a field of 2500 G is insufficient to fully decouple the relaxation. This indicates 
that at low temperatures there is an increasing electronic contribution, i.e. that as the 
temperature is decreased, magnetic fluctuations are slowing and falling into the muon time 
window. 

We have therefore fitted the zero-field muon polarization decay for different temperatures 
(examples are shown in the inset of Fig. [6]) with two components: a temperature indepen- 
dent Kubo-Toyabe function to represent the nuclear component, combined with a stretched 
exponential relaxation of the type exp((— At)' 3 ), where is a constant and A is the muon 
depolarization rate. The temperature dependence of A is shown in Fig. [6} It can clearly be 
seen that the depolarization rate begins to increase at T ^ 10 K, as the inverse susceptibil- 
ity becomes non-linear and the peak appears in the diffuse neutron scattering. There is a 
continual and relatively rapid increase until it reaches a plateau at T w 0.4 K, close to the 
proposed spin glass transition. 

D. Inelastic neutron scattering 

A limited inelastic neutron scattering investigation was also made. The results are shown 
in Fig. [7} The peak in the magnetic scattering previously described and attributed to 
nearest neighbor AFM correlations is at \Q\ m 1.25 A^ 1 . In Fig. [7] it can be seen that at 1.7 
K inelastic scattering extends from this position. There is no sign of any sharp feature, often 
visible in powder samples of ordered magnets, even when dispersion surfaces are averaged. 
The intensity is confined to this single feature, which is rather broad and weakly dependent 
on the energy transfer. Constant energy cuts further highlight the fact that the magnetic 
part of the inelastic signal can be fitted by the antiferromagnetic dimer structure factor at 
all energy transfers (Fig. [TJd). A constant |Q|-cut shows no indication of a gap (Fig. [7^). 
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FIG. 6: The muon depolarization rate in zero field as a function of temperature obtained from fits 
to a temperature independent nuclear response and stretched exponential electronic contribution. 
Typical fitted muon responses at different temperatures are shown in the inset. 

IV. DISCUSSION 

A spin liquid may be denned by analogy with an actual liquid, as a spin system with 
short ranged, gapless, dynamic correlations®!. We seek to make the case that despite strong 
quenched disorder and the previous description as a spin glass, LuCuGa0 4 is actually a spin 
liquid. We first discuss the structural disorder and secondly the liquid-like spin dynamics. 

In all the models fitted to the powder diffraction data the Cu 2+ and Ga 3+ ions were 
treated as randomly distributed on the triangular bilayers, as in all previous works on related 
materials. However, although the ordering of two different cation species on the triangular 
bilayer is frustrated (see Fig. [TJ), a completely random distribution seems unlikely: triangles 
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FIG. 7: a: S(Q,oj) at 1.7 K with Ei = 17 meV. b: Constant energy cuts through the feature at 
\Q\ rj 1.2 A" 1 visible in a, with fits to the nearest neighbor AFM structure factor. The peak around 
Q = 1.2 A" 1 fits well to this model at all energies. There is an increasing phonon contribution at 
higher \Q\. c: Cut at \Q\ = 1.2 A" 1 and 1.9 A" 1 (background). 
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containing entirely Cu 2+ or Ga 3+ will incur a Coulombic cost. We therefore advance the 
suggestion that in these materials the cations must be correlated. The nearest neighbor part 
of the Coulomb energy will be minimized by configurations of the two cations which map to 
the groundstates of an Ising antiferromagnet on the same lattice, as in other systems with 
two cation species on a frustrated latticed Structural diffuse scattering would be expected as 
a result of these short range correlations. Unfortunately, the scattering lengths of copper and 
gallium are extremely similar (bc n = 7.718 barn and &Ga — 7.288 barn) so direct detection 
of short range cation order is effectively impossible. 

The importance of this point lies in the resulting quenched disorder configuration. In 
the case of completely random cations, the topology of the spin system can be compared 
to percolating clusters, as in Ref. [181 In the case of correlated cations, it appears that the 
topology of the spin system would take the form of branching loops^. Less direct evidence 
exists for the formation of short-range correlations amongst the cations. For example, one 
can argue that the existence of two defined positions for the Lu 3+ ions is due to the mini- 
mization of electrostatic interactions amongst both the layer and interlayer cations. Without 
direct evidence for cation correlations, we conclude that the copper and gallium are certainly 
disordered, and that there is a possibility that they realize a non-trivial network structure. 

We turn now to the magnetic correlations. Although we have fitted the peak in the 
polarized neutron scattering data with the so-called dimer structure factor, we note that 
this is actually a general description of purely near-neighbor antiferromagnetic correlations. 
The inelastic scattering data shows that at 1.5 K dynamical correlations exist. The dif- 
fuse scattering observed in the polarized neutron scattering experiments may be regarded 
as integrating over the quasi-elastic part of these fluctuations. The temperature indepen- 
dence of the diffuse peak implies that the correlations remain dynamic at all temperatures 
measured: if the spins became significantly more static, one would expect to see transfer 
of spectral weight into the elastic line, certainly resulting in stronger diffuse scattering and 
possibly a sharper peak shape. The literal interpretation of the dimer structure factor is 
further precluded as there is no sign of significant singlet formation: the susceptibility does 
not decrease, there is no clear sign of a singlet-triplet gap in the inelastic neutron scattering 
data, and the //SR data implies that the local field in the sample is finite and fluctuating at 
all measured temperatures. 

Canonical spin glasses such as AgMn(0.5 at. %) have been extensively investigated by 
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/iSR 32 . The depolarization rate increases as the system approaches T g but then falls rapidly. 
In contrast, in LuCuGa04 the depolarization rate increases steadily as the spin correlations 
build up, but there is no evidence of freezing. We therefore also rule out the existence of a 
spin glass transition. 

The literature provides various examples of similar neutron scattering responses, for ex- 
ample energy independent continua at a particular \Q\ can be found in deuteronium jarosite 
((D 3 0)Fe 3 (S04)2(OD) 6 , kagome lattice) 30 ; SCGO (SrCr 8 _ x Ga4 +x 19 , x = 0.87, pyrochlore 
bilayer)P^ and Yo.5Ca .5BaCo 4 07 (kagome)^. However, the behaviour of LuCuGa0 4 is 
most strikingly similar to that of ZnCu3(OD) 6 Cl2 (Herbertsmithite), which is thought to be 
a structurally perfect s — 1/2 kagome antiferromagnet. In Herbertsmithite, x _1 curves be- 
low the Curie-Weiss straight line, which may be attributable to Zn-Cu disorder^ESl. Detailed 
/iSR investigations^ found a stretched exponential relaxation and increasing depolarization 
rate at low temperatures which was attributed to slowly relaxing defects, with the remain- 
der of the spins fluctuating too rapidly to be observed by /zSR. The neutron scattering 
response^ consists of a non-dispersive, gapless feature extending across all observed ener- 
gies at the wavevector transfer characteristic of the near neighbor distance. It can be fitted 
by the dimer structure factor at all available energies and measured temperatures between 
2 and 120 K. Detailed investigation of the behavior of x"(QjT) extracted from the neutron 
scattering data shows that it does not have the expected E/T scaling expected in theories 
of spin liquids (this is also the case in deuteronium jarosite®'). While we do not have suffi- 
cient data to compare with these observations, we suggest that the overall similarity of the 
two compounds means that the same physics is realized and may be a general property of 
two-dimensional, s = 1/2 spin liquids, which is robust to the presence of disorder. 

V. CONCLUSION 

We have extended the characterization of LuCuGa04 and shown that it is not a spin glass, 
but a spin liquid. Although a model of randomly disordered cations is sufficient to explain 
diffraction data, this may be due to the almost identical scattering lengths of copper and 
gallium and further investigation of the topological properties of these materials promises 
to be interesting. At moderate temperatures the spin correlations are shown to be dynamic 
and gapless. There is a complete absence of spin freezing down to the lowest available 
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temperature in LuCuGa04, suggesting that the dynamic correlations persist despite the 
structural disorder. LuCuGa0 4 therefore provides another opportunity to investigate the 
scaling behaviour of x(Q,T), demonstrated to be important in other two-dimensional spin 
liquids. 
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